Frequency-selective near-field radiative heat transfer between photonic-crystal slabs: 
A computational approach for arbitrary geometries and materials 
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We demonstrate the possibihty of achieving enhanced frequency-selective near-field radiative heat 
transfer between patterned (photonic crystal) slabs at designable frequencies and separations, ex- 
ploiting a general numerical approach for computing heat transfer in arbitrary geometries and mate- 
rials based on the finite- difference time-domain method. Our simulations reveal a tradeoff between 
selectivity and near- field enhancement as the slab-slab separation decreases, with the patterned 
heat transfer eventually reducing to the unpatterned result multiplied by a fill factor (described by 
a standard proximity approximation). We also find that heat transfer can be further enhanced at 
selective frequencies when the slabs are brought into a glide-symmetric configuration, a consequence 
of the degeneracies associated with the non-symmorphic symmetry group. 

PACS numbers: 



Radiative transfer of energy from a hot to a cold body 
is well known to be enhanced (even exceeding the black- 
body limit) when the bodies are brought close enough 
for evanescent fields to contribute flux [1-9]. In this 
paper, we demonstrate that near-field radiative trans- 
fer can be greatly modified by using periodically pat- 
terned photonic-crystal (PhC) surfaces, with frequency- 
selective enhancement that can be controlled by choosing 
the geometry (rather than relying on material or plas- 
mon resonances available only at long wavelengths 
[l2|). Until now, investigations of near-field transfer in 
microstructured geometries have been hampered by the 
lack of computational modelling techniques, and we em- 
ploy a new rigorous approach based on directly simu- 
lating Maxwell's equations in time with random thermal 
sources, an extension of the Langevin approach we previ- 
ously used to model the emissivity of a single body [l3| . 
Previously, aside from semi-analytical results for planar 
structures 111, Q, [HI , formulations have been developed 
that in principle handle arbitrary geometries but which 
thus far have only been evaluated using Fourier methods 
specialized only for pairs of spheres and/or plates [l6l - 
[ill . Multilayer planar structures are tractable [22|, but 
(Fabry-Perot) resonant modes created in such structures 
are not evanescent in the air gap, unlike leaky modes 
in transverse-patterned structures (each of which has 
both evanescent and propagating fields in the gap), so 
multilayer films do not lend themselves to frequency- 
selective near-field enhancement. In the far field, it is 
known that more complicated structures such as PhCs 
can be designed to resonantly enhance radiative trans- 
fer at desired frequencies 23-25|, which is crucial for 
applications such as thermophotovoltaic cells in which 
only certain frequencies can be efficiently converted to 
power [sl, [25l - [27j . One would like to obtain similar 
short- wavelength (< 2 /xm [25|) enhancement of near- 



field effects, whereas plasmon resonances are too long- 
wavelength for such applications. In a simple model sys- 
tem consisting of two PhC slabs, thin films with peri- 
odic grooves [Fig. [TJtop)], we show that the resonant 
leaky modes created by the periodicity yield orders of 
magnitude enhancement in the flux at designable wave- 
lengths even for moderate separations (100s of nm to 
microns for infrared wavelengths) compared to similar 
structures in the far field, starting with weakly absorb- 
ing thin slabs that transfer < 1/1000 of the flux between 
black bodies in the far field. Furthermore, we show that 
the selective enhancement can be almost doubled at se- 
lective peaks by using a glide-symmetric configuration 
[Fig. [Sl^bottom)], due to degeneracies resulting from the 
properties of the non-symmorphic symmetry group (28| . 
Ultimately, there is a tradeoff between enhancement and 
frequency selectivity — much larger enhancement is the- 
oretically attained for nm separations [4, 5, 7-9] where 
geometric resonances have no effect, at the cost of fre- 
quency selectivity and much more difficult fabrication — 
and this letter offers a first glimpse of the practical design 
space that is available to optimize these considerations. 

Given two arbitrary bodies, at temperatures Ti and 
T2, their radiative heat transfer i^(co', Ti, T2) is: 

H{co;Ti,T2) = $H [e(w,Ti) - e(w,T2)] , (1) 

where <l> is the flux into a single object due to random 
(white noise) current sources present in the other object, 
and O is the Planck distribution [il. (Physically, for 
linear electromagnetism it is equivalent to use Planck- 
distributed thermal fluctuations or to use white-noise 
sources multiplied afterward by O.) Computationally, 
this formulation can be expressed directly in terms of 
a Langevin model [HI, in which white- noise sources are 
introduced into the evolution of Maxwell's equations in 
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FIG. 1: Top: Schematic geometry of two PhC slabs of thick- 
ness h = 0.2a, separated by a distance with periodic air 
grooves of period a and width t = 0.2a. Bottom: Modal fre- 
quencies uj(kx) of the TE (H • z = 0), kz =0 modes of an 
isolated slab; note the presence of both leaky (gray region) 
and guided modes. Insets: mode profiles at kx — {0,7r/a}. 



a finite-difference time-domain (FDTD) method. (More 
sophisticated approaches for solving this problem exist, 
e.g. non-stochastic methods in the frequency domain, 
but the sacrifice in efficiency of our FDTD approach is 
compensated by its generality and simplicity.) Although 
we previously only applied this method to equilibrium sit- 
uations [13], the extension to nonequilibrium situations 
is straightforward because the statistical independence of 
random currents in different objects allows one to calcu- 
late the flux due to thermal sources for one object at a 
time. Reciprocity allows one to calculate the flux <l> from 
one body to the other and infer the flux in the other 
direction, a fact that is implicit in Eq. ([T]). 

We investigate the proof-of-concept structure shown 
in Fig. [TJtop), involving two PhC dielectric slabs (thin 
films) of thickness h = 0.2a, separated by a distance 
d in the x direction, with z-oriented air grooves of pe- 
riod a and width t = 0.2a in the x direction. The 
permittivity of the slabs is taken to be of the Drude 
form e{uj) = Soo — (j/uj{uj + ^7), where Soo = 12.5, 
a = 0.2533 (27rc/a)^, and 7 = 1.5915 (27rc/a), approx- 
imating a dispersionless dielectric with Kee ^ 12.5 and 
Ime ^ 1 at wavelengths comparable to a. The compu- 
tation of <l> for this structure involves introducing uncor- 
related white-noise sources into the damped polarization 
equation [13] at each position (pixel) of one of the bod- 
ies [red slab in Fig. [TJ^top)], and integrating the resulting 



fiux spectrum over a surface surrounding the other slab 
[blue slab in Fig. [IJtop)]. The resulting spectrum is en- 
semble averaged over many simulations (~ 60) to reduce 
the noise level in the spectrum. We employ a periodic 
unit cell in x with Bloch-periodic boundaries (phase dif- 
ference e*^^'^), truncate the cell in the y direction using 
standard PML boundary conditions, and integrate <l> over 
The translation-invariant z direction can similarly be 
handled by integrating 2d simulations over kz- We first 
consider the purely 2d problem corresponding to kz = 
and comment on the full 3d problem further below. At 
the moderately large separations ~ a studied here (as- 
suming a in the microns), the TE modes (H • z = 0) of 
the structure exhibit stronger confinement than the TM 
modes [29], and end up dominating the heat transfer at 
all relevant cj, and therefore we only compute the TE 
contribution. 

Figure [T] plots the (TE, kz = 0) modal frequencies 
uj{kx) of a single (isolated) slab. Modes that lie in the 
gray region (above the u = ckx light line) are leaky 
modes [29] , which radiate into the air and therefore con- 
tribute to heat transfer between the slabs in the far field. 
Modes that lie below the light line are guided modes [loj , 
which evanescently decay in air and thus transfer energy 
only in the near field. 

Figure [2jtop) plots the flux enhancement r] = 
*^t,d/^o,oo as a function of a;, where ^t,d is the flux spec- 
trum between PhC slabs of groove width t and separated 
by distance (i, here normalized by the flux spectrum <l>o,oo 
of the unpatterned slabs {t = 0) in the far field {d 00) 
to clearly illustrate the effect of finite t and d. {^0,00 
is computed by Kirchoff's law from the reflectivity in 
the far field by transfer-matrix methods.) The <l>o,oo 
spectrum by itself is nearly flat, so the spectral features 
in Fig. [2] are not due to the normalization. Note that 
we are starting with very weakly absorbing slabs (thin 
films), i.e. ^0,00 is more than 1000 times smaller than 
the corresponding ^ = 1 of two black bodies. 

Radiative heat transfer can be significantly changed 
by the periodicity of the slabs [2^, as illustrated by 
the green curve in Fig. [2ftop), which shows rj for PhCs 
of groove width t = 0.2a in the far field. <l>o.2a,oo ex- 
hibits wide-bandwidth peaks, orders of magnitude larger 
than ^0,00- Not surprisingly, the bandwidths of these 
peaks match the bandwidths of the leaky-mode bands in 
Fig.lljbottom), whose integrated contribution over all ky 
leads to a smearing of the sharp spectral peaks that oc- 
cur at the leaky- mode frequencies for each ky. The sharp 
drop in r] between the peaks is a consequence of the pres- 
ence of pseudogaps between the leaky- mode bands. 

As the separation decreases, a number of dramatic fea- 
tures are observed in Fig.[2jtop), which shows 77 at three 
additional separations {d = 0.2a, d = 0.5a and d = a). 
At d = a, one can observe a dramatic increase in rj at 
low frequencies. This arises due to the evanescent cou- 
pling and contribution of guided modes to h. This near- 



3 



10^ 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 
frequency co (27ic/ a) 



0.16 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 



frequency CO {2nc/a) 

FIG. 2: Top: Flux enhancement 77 = $o.2a,d/^o,oo as a func- 
tion of cj (units of 27rc/a), where ^t,d is the flux spectrum 
between PhC slabs of groove width t and separation at 
various d = 0.2a, 0.5a, a, oo. Solid (dashed) lines denote 77 for 
PhC (unpatterned) slabs. Bottom: Flux spectrum $0. 2a, 0.2a 
[units of hc{27v/a)'^] for a symmetric (s = 0, red line) or glide- 
symmetric (s = 0.5a, black line) slab-slab configuration. 



field coupling increases further at smaller and addi- 
tionally (especially in the d = 0.2a curve) one can see 
a dip corresponding to the gap in the guided modes of 
Fig. [D^bottom), and sharp peaks at the edges of the gap 
(a consequence of van Hove singularities [29|). There is 
also an enhancement of <l> at small d at the leaky-mode 
peaks because, in a patterned slab, leaky- modes have an 
evanescent as well as a radiative component [29| , and the 
former enhances the coupling at short distances. As for 
the guided modes, this evanescent coupling is further en- 
hanced by the van Hove singularities of the zero-slope 
"slow-light" band edges. Also, the leaky modes split 
into bonding and antibonding pairs as d decreases, and 
a similar splitting is visible in the guided mode peaks at 
even smaller a phenomenon visible even more clearly 
in Fig. [2jbottom). 

A striking phenomenon in the near-field interaction oc- 
curs when one patterned slab is shifted in the x direc- 



tion by 0.5a with respect to the other slab, as shown in 
Fig. EJbottom). In this case (plotted here as absolute <l> 
on a linear scale), the bonding-antibonding peaks of the 
symmetric configuration (solid red curve) merge at cer- 
tain cj, and can even add to approximately double the 
heat transfer of the unpatterned slab. This arises due to 
the glide symmetry of the shifted structure, whose non- 
symmorphic symmetry group supports degeneracies at 
the Brillouin zone edge {kx = n/a) [2^, causing the cor- 
responding bonding-antibonding mode splittings to dis- 
appear. In contrast, the frequency splitting correspond- 
ing to the kx = band edges, e.g. at ~ 0.9 (27rc/a), 
persists even for the glide-symmetric structure. Simi- 
larly, the gap in the guided modes of Fig. [TJbottom), 
near ~ 0.4 (27rc/a), is greatly modified by the presence 
of the other slab in the mirror-symmetric case (red curve) 
but not in the glide-symmetric case (black curve). 

Figure [2](top) also shows the enhancement for unpat- 
terned {t = 0) slabs. These also exhibit a striking near- 
field enhancement for d < a, but in this case the enhance- 
ment is broad-band, because there are no standing-wave 
leaky-mode solutions in an unpatterned slab to induce 
frequency selectivity. In particular, at relatively large 
separations (e.g. d = a), the patterned slab exhibits more 
than 10 times the enhancement of the unpatterned slab, 
due to these resonant contributions. At d = 0.5a, the un- 
patterned slab is better for small uj and is worse for larger 
cj, where the latter is due to the fact that the evanescent 
interactions become shorter-range at high uj and vanish 
in the absence of resonant enhancement. At d = 0.2a, 
the unpatterned slab has larger enhancement over the 
whole bandwidth shown here, but lacks frequency selec- 
tivity. In general, as d decreases, the periodic structure of 
the PhC slabs ceases to matter, and the enhancement of 
the peaks relative to the overall near-field enhancement 
is lost. At even smaller not shown here, the interac- 
tion between the PhCs can be described by a proximity 
approximation in which adjacent surfaces are treated as 
parallel plates 0, 0i| multiplied by a fill factor. 

Figure [3] plots the total (integrated) heat transfer Ht^d 
as a function of separation d for both patterned (i^t,d, 
squares) and unpatterned (i^o,^, thin lines) slabs of pe- 
riod a = 2/im, at various temperatures T2, and Ti = 0. 
At T2 = 500 K, dominated by low-cj contributions, 
the unpatterned slabs yield larger H over almost ev- 
ery d (even as d 00), which increases for small d 
due to near-field enhancement. Moreover, we find that 
^ = Ht^d/Ho^d 1 — t/a = 0.8 as d ^ 0, as expected 
from the proximity approximation. (For TE modes, H 
asymptotes to a finite value as d ^ 0, whereas for TM 
modes, whose contributions are negligible at the separa- 
tions we consider here, H eventually diverges At 
T2 = 1000 K, larger uj begin to dominate, and the pres- 
ence of PhC peaks [Fig. [2] causes a significant increase 
in H. In the (purely radiative) limit d ^ 00, the PhC 
slabs perform much better than the unpatterned slabs 
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~ 3), although the tradeoff between enhancement and 
selectivity is apparent as there exist a ~ 7 fim be- 
low which ^ < 1 (eventually 0.8). As expected, dc 
decreases as T2 increases due to the increasing role of 
higher- frequency resonances. At T2 = 2000 K, we find 
that ^ ^ 5 at larger increases at intermediate d and 
falls below 1 at ^ 2 /im. 

Figure [3] illustrates the advantages of using PhCs at 
large and intermediate d^ where both selectivity and near- 
field effects coexist. However, in many applications, the 
figure of merit is not the overall H but the integrated 
H over a finite frequency bandwidth. In thermophoto- 
voltaics, heat transferred from a hot body to a cold 
semiconductor is converted to electricity via electronic 
bandgap transitions, which means that photons of fre- 
quency smaller than the bandgap frequency ujg (bandgap 
wavelength of a few microns) are lost in the conversion 
process [9[. The left inset of Fig. [3] shows a correspond- 
ing ad hoc figure of merit: the windowed heat trans- 
fer Hw = / duj H{uj^Ti^T2)n{uj) between slabs of period 
a = 1/im, as a function of separation and at tem- 
peratures Ti = 300 K and T2 = 1500 K. The window 
function n is similar to the "quantum efficiency" func- 
tion of a typical semiconductor of bandgap wavelength 
27Tc/ujg = 2.5 /im [9]. The spectra H{uj) and n{uj) are 
plotted on the right inset of Fig. [3] at d = 2 jam. Notice- 
ably, the presence of high-frequency peaks in the PhC 
makes its over an order of magnitude larger than 
that of the unpatterned slab for d > 1 /im; the situa- 
tion is even more favorable for the PhC because the low- 
frequency photons in the unpatterned case are wasted. 

Our calculations of H were restricted to kz = for 
computational convenience. However, preliminary cal- 
culations of the /c^-integrated H{uj)^ for a single separa- 
tion, reveal no qualitative changes in the spectrum, ex- 
cept for a slight broadening of the spectral peaks. As 
noted above, the geometry of Fig. [T] is in no way optimal 
and represents only a proof of concept: we believe that a 
2d-periodic geometry will yield even better performance, 
an exciting direction for future numerical studies. 
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FIG. 3: Heat transfer Ht,d [units of /ic(27r/a)^] between un- 
patterned (t = 0, lines) or PhC (t = 0.2a, squares) slabs of 
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Ti = 300 K and various T2. Left inset: hw as a function of 
d (units of a) for slabs of period a = l/xm, and T2 = 1500 K, 
where is computed by windowing the spectrum h{uj) 
with a function that excludes contributions of wavelengths 
> 2.5 /xm. Right inset: Ht^d (gray) and ifo,d (black) at a 
typical separation d — 2 /xm, and the window function n{(jj) 
(orange line). 
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